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We argue that the curvaton decay takes place most naturally by way of a broad parametric resonance. 
The mechanism is analogous to resonant inflaton decay but does not require any tuning of the 
curvaton coupling strength to other scalar fields. For low scale inflation and a correspondingly low 

- - - mass scale for the curvaton, we speculate on observable consequences including the possibility of 

QP \ stochastic gravitational waves 

o 

^ ^ , INTRODUCTION 

D ■ 

^y^ . The curvaton scenario [1-4] provides an alternative mechanism for producing the primordial density perturbations. 

The curvaton is a Hght field which is subdominant during inflation but whose energy density will eventually become an 
important (if not actually the dominant) energy component in the universe. Because of its lightness, during inflation 
the curvaton acquires fluctuations which at the time of curvaton decay are transformed into curvature perturbations. 
So far, the decay of the curvaton has mostly been treated as a perturbative process, see e.g. [5, 6], where the 
curvaton condensate decays into relativistic particles, determined by a rate F and occurring roughly when H ^ T. 
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I . This happens at the stage when the curvaton field is oscillating about the minimun of its potential, whence it behaves 
O ' like pressureless dust so that its relative energy density grows with respect to radiation produced by the inflaton 
-*— > I decay. The curvature pertubation induced by the curvaton decay is mostly adiabatic but can give rise to a large 
jrt non-gaussianity [6, 7] if the curvaton energy density is too small^; this implies a lower limit on the fraction of the 
I— '. energy density in the curvaton condensate at the decay time. Another generic feature of curvaton models with no 
^ I mixed inflaton and curvaton perturbations is the small amplitude tensor perturbations [2, 6]; this is a consequence of 
<^ . the low inflationary scale required to keep the inflaton perturbations negligible. ^ 

(S) ■ However, since the curvaton is a coherently oscillating condensate, it is possible to excite other bose flelds in 

00 . a regime of parametric resonance, in a manner similar to inflationary preheating [10, 11]. As a consequence, the 
^J^ ' curvaton decay process could have important cosmological ramiflcations for e.g. non-gaussianities, which are known 
^— ; . to be enhanced [12, 13] in some preheating scenarios, as well as for stochastic gravitational waves, as will be discussed 
t^^ ' below. Preheating is a violent process in which the occupation numbers of flelds coupled to the curvaton grow 
^^ , exponentially. In effect, dense lumps of decay products are created very rapidly and randomly in the real space, 
and as these highly inhomogeneous concentrates of energy density collide, there arise large-amplitude, high-frequency 
gravitational waves, as has been discussed recently [14, 15]. These do not leave a signal on the CMB but could possibly 
be observed in current or future gravitational wave detectors such as LIGO [16], LISA [17] or BBO [18]. 

In the present paper we study the possibility of preheating in generic curvaton models and explore the main 
$H ' consequences of such a scenario (preheating has been discussed in the context of MSSM flat direction curvatons in 
[19]). Indeed, the curvaton has to couple to other flelds if it is to decay; if it couples to other scalars, then there 
would be two options. Either the additional scalars would be light and shifted from their minima during inflation, in 
which case they would simply be additional curvatons, or they could be massive and driven to their minima. In the 
flrst case one still would have to flnd a way for the curvaton(s) to decay. However, in the latter case there can arise a 
parametric resonance during curvaton oscillations with the ensuing rapid decay of the curvaton. 

We flnd that the preheating process in the curvaton scenario is qualitatively similar to the preheating after inflation 
although the resonance dynamics can be somewhat altered. This is because the curvaton, unlike the inflaton, could 
remain subdominant during the entire preheating stage. Moreover, in the scenario studied here the backreaction 
of the excited bose flelds always becomes signiflcant at the end of preheating whereas in the standard inflationary 



00 
O 



X 



a 



^ strictly speaking, this holds only for a quadratic curvaton potential, see [8]. 

^ See however [9] where it is demonstrated that the amplitude of CMB-scale tensor perturbations in the curvaton model can be directly 
proportional to the amount of the produced non-gaussianity due to second order effects. 



case the preheating may terminate already at an earlier stage. Since the gaussian part of curvaton perturbations on 
superhorizon scales is not affected by preheating and the final decay of the curvaton will eventually take place in a 
perturbative reheating process, the possible preheating stage does not alter the standard predictions of the curvaton 
model [6]. However, interestingly enough we find that preheating may generate stochastic gravitational waves with an 
amplitude much larger than the infiationary tensor perturbations in the curvaton model, and we discuss the feasibility 
of observing such signatures of curvaton preheating. The study of the non-gaussianities is left for a future work. 

THE INFLATIONARY PERIOD 

Preheating, or resonant decay of the curvaton, may occur in a large class of curvaton models with several interacting 
scalar fields, as we will now discuss. To illustrate the main features of the process, let us consider a very simple model 
with the potential 

V = V{^) + V{a,x), (1) 

where (j) is the infiaton, a is the curvaton and x is an additional scalar field which the curvaton decays into. Note that 
X should not be yet another curvaton. We assume that the infiationary dynamics is completely determined by the 
infiaton field, H^, = 7?* (</>), which requires that the energy densities of the curvaton and the x field are much smaller 
than H^Mp. We also assume that all fields have canonical kinetic terms. 



Let us further assume that 



na,x)-im^.2 + igVx^ + lAx\ (2) 



where 5 < 1 and A < 1 are some coupHng constants. For certain field and coupHng values, both a and x could act 
as the curvaton but this is not the case we wish to study here. We therefore concentrate on the regime where the 
effective mass of the x field is larger than the Hubble parameter ga > iJ, and the contribution to the curvaton mass 
coming from the x field can be neglected x ^ ma/ g. In this case the field x will not acquire perturbations during 
infiation and it simply represents a generic bose field that is coupled to curvaton but will not itself act as a curvaton. 
Under these assumptions, the value of the x field is driven towards zero exponentially fast during infiation. 

The curvaton is effectively massless during infiation mo- <C i?* and acquires perturbations proportional to the 
Hubble parameter 5a ~ iJ*. In the standard scenario where the curvaton is assumed to decay perturbatively at 
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some stage after infiation, the amplitude of primordial curvature perturbations can be estimated by P^ ~ rH^,/a 
[2, 6], where r ^ Pa/ P measures the fraction of the curvaton of the total energy density at the decay time. Since the 
curvaton is very light, its equation of motion is approximatively given by 3-ff»CT ~ —m^a (for a discussion of curvaton 
dynamics in various potentials, see [20]). If infiation lasts long enough, the curvaton will eventually roll down to the 
region in its potential where the quantum effects overcome the classical drift, \d\ < H^ . The typical curvaton value 
during infiation in this case is thus given by 

2 



H*[?^] ■ (3) 



In the Hmit r ^ 1, this yields an estimate 
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(4) 



for the curvaton contribution to the total energy density at the beginning of oscillations. 

THE PARAMETRIC RESONANCE 

Preheating in the curvaton model can take place when the curvaton starts to oscillate after infiation and gives a 
time dependent effective mass for the x field that initially has a vanishing vev. The situation is analogous to the 
standard case of preheating in the conventional infiationary models [10, 11]. As usual, an efficient preheating requires 
broad resonance bands which implies the condition 



at the beginning of oscillations. In our case, this necessarily holds as a consequence of the mass condition for the x 
field, g > tt^ja, combined with the lightness of the curvaton mc <C -ff*. 

Depending on the fraction of the energy density in the curvaton field at the time of oscillations, denoted as rose, 
there are two different possibilities for realizing the preheating stage. If the curvaton dominates the energy density 
at the beginning of oscillations with rose ~ 1, preheating proceeds exactly in the same manner as in the standard 
infiationary case. In the opposite limit, rose ^ 1, the universe is not dominated by the oscillating curvaton but by 
the infiaton decay products assumed to be homogeneous radiation. This however leads only to minor differences as 
compared to the standard case [11] as we will discuss below. 

In the Hmit rose ~ 1 the universe is effectively matter dominated during the curvaton oscillations. The time 
evolution of the curvaton is given by 

cr(t) — Sin(77lcrt) = G svi\(m„t) , (6) 

where a ex t~^ measures the amplitude of oscillations. For definiteness, we choose to = T^l^^fn^) as the beginning of 
oscillations and set ao = 1. The initial amplitude is then given by a = 2(7*/7r, where ct* is the value during infiation. 
Using a rescaled field variable X\^ — a?''^Xk, the equation of motion for the x field in Fourier space reads 

Xk+(^+9^^^\xk^0, (7) 

where terms proportional to H^ <C g^<J^ and the small contribution from the quartic term in (2) have been neglected. 
The solutions can be formally written as 

X,(t) = ^e-r-'^* + M£e^r-^* , (8) 

\/2lo \/2uj 

and the occupation number of a given mode is determined by n^ = |/3fcp. In the broad resonance regime, the 
occupation numbers rik remain constant for most of the time since the mass of the x field is slowly varying and much 
larger than the curvaton mass nif, [11]. However, as the oscillating curvaton crosses the origin during each cycle of 
oscillation, the x field becomes effectively massless and its frequency Wfe — ^Jikja)"^ + g'^a'^ changes in a non-adiabatic 
fashion ijjk ^ w^. This can lead to copious production of x particles via a resonant decay of the curvaton condensate. 
The particle production takes place at very short time intervals around times i„ defined by a{tn) = 0. The curvaton 
value in the vicinity of these intervals can be approximated by \a\ ~ anfna{t — tn)- By substituting this into (7), the 
equation of motion becomes 

^ + {>,l + rl)Xk = Q, (9) 



where we have introduced the variables /c* = y^gm^a, Kn — k/{ak^), r„ = k^{t — tn) and all the time arguments not 
expHcitly written out are understood to be evaluated at t„. The solutions for this equation are paraboHc cyHnder 
functions W^(— k^/2; ±-\/2t„). By comparing them in the subsequent adiabatic regimes i„_i < f < t„ and t„ < t < 
tn+i one can estimate the change of occupation numbers n^ around the times tn when the system is non-adiabatic. 
In the limit of large occupation numbers, n^ ^ 1, the result is given by [11] 

<+i = < exp(27rM^) , (10) 



where 27r^^! = ln(l + 2e^^'^^ — 2sin0„ e^^'^"/^v 1 + e^^*^" ) and On is a phase factor that is essentially a random 
variable for each resonance time tn- This is why preheating in an expanding space has a stochastic nature. 

The system formally describes scattering of waves in a parabolic potential. Each resonance time t„, when the 
curvaton hits the origin, corresponds to a scattering event and changes the occupation numbers of the x modes. The 
changes can be significant in the regime ttk^ < 1 which corresponds to the momentum band 

^^ < kl{tn) ^ gm^an ,^^. 



On average, the occupation numbers grow exponentially in this band, although the exponent ^fc in (10) can also be 
negative for some scattering events f„, and the comoving particle number density as a function of time can estimated 
as 

n{t) = -—--; / d^k ndt) - 10'-'^ , * =6^^"-* . (12) 



Here /i is the maximum of the effective growth index /i|,*^ obtained from equation (10) by defining ii'^t = J dt/ife(i). 
As the modes of the x field grow exponentially, the term g^x^cP' in the curvaton potential becomes comparable to the 
mass term m^a^ and the backreaction effects which eventually terminate the resonance need to be taken into account. 
The backreaction becomes important [11] when the number density of x particles (12) grows to n(ibr) ^ "^a^ 1 9 which 
yields an equation 



*br - T^ln ( i^iip^ ) , (13) 



from which tbr can be solved. The condition for broad resonance (5) becomes violated when gcr{t{) ~ iria-, which 
yields 

^-^>^. (14) 

The last inequality follows from the condition g > H/a which is required to keep x massive during infiation. The 
backreaction generically takes place much before the end of broad resonance since the above equations yield ibr ^ if 
for nia/H^ < 10~^, assuming the curvaton value during inflation is given by equation (3). The interaction between 
the curvaton condensate and the x particles breaks the condensate into curvaton particles which eventually come to 
dominate over the x particles [11]. As usual, preheating needs to be followed by a reheating stage which eventually 
yields the hot big bang era. The simplest assumption [11] is that preheating just prepares slightly different initial 
conditions for the reheating stage during which the curvaton decays into standard model flelds. In this case the 
standard predictions of the curvaton scenario for the gaussian part of density fluctuations remain largely unaltered. It 
is possible though that the non-thermal epoch after preheating could signiflcantly affect the picture and lead to serious 
constraints on the curvaton scenario; this question should be addressed numerically by studying the thermalization 
process. However, even if the gaussian density fluctuations would not be affected, curvaton preheating leads to 
production of gravitational waves which could be detectable as discussed below. In addition, although we do not 
discuss this possibility further here, we note that preheating could also generate a signiflcant level of non-gaussian 
density perturbations as can happen in the usual inflationary case [12, 13]. 

If the curvaton is subdominant at the beginning of oscillations, rose ^ 1, the above discussion gets slightly modified 
but the qualitative results remain the same. The exact time evolution of the curvaton in the radiation dominated 
background is given by Bessel functions but the solution is well approximated by 

"(^)" (m.^ + f)3/4 ^^"("^-^+i) ■ (15) 

This differs from the matter dominated case (6) only by a phase and the slightly slower decrease of the amplitude 
a ex t~'^/^. It is therefore rather evident that the analysis of preheating will not significantly differ for the matter 
dominated case. The equation of motion for the x field can be expressed exactly in the same form as equation (9) 
except that the time evolution of the variables k„ and r„ differs from the matter dominated case due to the different 
evolution of the scale factor a{t) and the curvaton amplitude a. This does not affect the stochastic nature of the 
successive scatterings at times t„ and the occupation numbers of the x field grow essentially at the same rate as 
before. We therefore expect that the number density of x particles in the rose *C 1 Hmit can be estimated by the 
same equation (12) as in the case rose ~ 1. This conclusion can be checked by solving the equations of motion for the 
curvaton and the x field numerically. As illustrated in Fig. (1), the results in the matter and radiation dominated 
cases indeed coincide quite precisely in the regime where the backreaction of x particles can be neglected. The time 
when the number density of the x particles has grown large enough to yield a significant backreaction can be estimated 
as before using equation (12) and the condition n(ibr) ~ fn^^/g. Due to the smaller suppression by the scale factor 
in (12), the backreaction will occur slightly earher than in the matter dominated case, as also seen in Fig. (1), but 
the effect is very small. In analyzing the dynamics after the onset of backreaction, the expansion of the universe can 
be neglected [11] and the radiation domination should not affect the situation here either. 

GRAVITY WAVES 

Since the decay of the curvaton can take place in a regime of parametric resonance, as was demonstrated above, it 
is natural to expect that gravitational waves will also be produced when it occurs. Indeed, the physics is similar as 
in the infiationary preheating: an oscillating field excites the production of particles in other bose fields at the times 
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FIG. 1: Comparison of resonant decay of the curvaton in radiation and matter dominated backgrounds. The figure shows Inn 
for the X field in the cases r=l (dashed lines) as well as r < 1 (solid lines). It also shows the limits above which backreaction 
becomes important in both cases. The parameters for these plots are (Mp = 1), m = 10~^°, g = 10~*, and a* — 0.0002. The 
onset of backreaction takes place after about 30 - 35 oscillations of the curvaton which corresponds to growth index /i ~ 0.07. 



when their effective frequencies change non-adiabatically. According to [15], in such circumstances the generation of 
gravitational waves is inevitable. In what follows we make extensive use of the analysis in [15] which can be directly 
applied to our case with appropriate modifications. 

Let us first reiterate the salient features of production of stochastic gravitational waves, which can be represented 
as the transverse-traceless part of the metric perturbations 



ds'^ = a'^{ri) [-drj'^ + (% + h,j) dx'-dx^] 



(16) 



where dihij = h 
momentum tensor T,^-'" ^ 



and di] = a ^dt. Gravitational waves are sourced by the transverse-traceless part of the energy 
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d^^h,{k) +ik'-^\ hj{k) = 167rGa(7y)r^^(k) , 



(17) 



where hij — ahij. In the case of radiation domination 9^a = so equation (17) can be solved with the appropriate 
Green functions of the harmonic oscillator. Since we will be considering modes with k :^ aH the same applies to a 
good approximation for the matter dominated case. The energy density in gravitational waves for subhorizon modes 
can then be written as [15] 
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Considering the gravitational waves being produced up to a time ?]i we have 
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where Sk for a universe filled with a scalar field x reads 
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^ We are considering here classical gravitational waves generated by the motion of matter and not vacuum fluctuations amplified by the 
cosmic expansion. 



Here (k, p) is the angle between the unit vectors and the factors sin(fc77) and cos{kr]) come from the harmonic oscillator's 
Green functions. Taking into account the expansion of the universe from rji, the present day energy density spectrum 
of the gravitational waves becomes 

n,^h^ = 4^M (ll.) ^ 0,0/.^ , (21) 

where we have assumed radiation domination, rose ^ 1, throughout the generation of gravitational waves for simplicity, 
g refers to the number of relativistic species and flroh^ is the energy density of radiation today. 

According to equation (11) the typical physical momenta amplified during preheating are /c* ^ m^q^'^ and the 
gravitational wave signal will also be peaked around A:,. The corresponding frequency today is [15] 

A ^ 4^ 10^° Hz ~ g'/'rli^af/' 10^° Hz (22) 

where the subscript "p" refers to the time when the peak was formed and w is the equation of state parameter, 
w = 1/3 for rose ^ 1 and w = for rose ~ 1- The scale factor at the beginning of oscillations has been normaHzed to 
unity as before. Using equation (4), which holds for models where the curvaton value during inflation is given by (3), 
and taking into account the constraint (/ct* 3> i?* required to keep x massive, one finds the bound 



The frequencies accessible to direct detection experiments range from / ~ 10"'* Hz (LISA) to / ~ 10"^ Hz (LIGO). 
According to equation (23) , the gravitational waves produced in the curvaton preheating can fit this observable window 
only if the infiationary scale is low enough. If the curvaton would decay perturbatively and if there are no significant 
changes in its effective mass after inflation, the inflationary scale is constrained by iJ* > lO^GeV [21]. This would yield 
Ik ^ lO'^Hz which unfortunately sHghtly escapes the observable window. However, equation (23) is a rather crude 
estimate and a proper numerical computation is required to really flnd out the actual frequency of the gravitational 
waves. Moreover, the bound i?* > lO^GeV does not directly apply here since the effective curvaton mass increases 
during the non-linear stages of preheating and the highly non-thermal epoch after the end of preheating may also 
significantly affect the expansion history of the universe as compared to the perturbative curvaton decay. 

Having discussed the frequencies of the gravitational waves, we are now ready to return to the question of their 
amplitude. Concentrating first on the linear regime where the backreaction of x particles can be neglected, we may 
use equation (20). Gravity waves will only be produced when the frequency of the x modes varies non-adiabatically, 
near the zeros of the oscillatory curvaton. Near those zeros the equation of motion for the x field is given by (9) 
(remember X = a^'^x) and its solutions are parabolic cylinder functions W{—k'^/2] ±v2r„) 

X," = AIW{^kI/2- V2t„) + i?IW(-4/2; -V2r„) . (24) 

This may be matched to the adiabatic solutions (8) in the region where the two solutions overlap. Assuming that the 
occupation numbers (10) always grow with the fastest possible rate, this yields up to a phase factor [15] 
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W — -^VK(-..y2;-V2r„), (25) 

where ^^ = Vl + e^^*^^ — e^'^'^ '^. Using this form and switching to the comoving time in (20) we find 

1 
Sk{n) ^ ^^^'^^^^^ K^ j du{l- u'f I dPP'^p gi^.pi Up (Tf )n|i,,p| (Tf ) (/2 + /2) , (26) 

-1 

where K = k/iria- and Is, /c are integrals of the parabolic cylinder functions. Except for the prefactors in front of the 
integral and the exact behaviour of the occupation numbers, this result coincides with equation (78) in [15] derived 
for the gravitational waves produced in the Hnear regime of preheating after X(j)'^ inflation. 

Ignoring the differences in Uk for a while, we may then relate the amplitude of the gravitational waves (21) obtained 
in [15], (r^gw/i^)^ -4, to the amplitude in our case as 

f^gw/.^-10^^(^)'(^gw/.^).,. , (27) 
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where subleading contributions coming from different values (and powers) of the scale factor and q have not been 
included. For a curvaton mass at the TeV scale, the amplitude in (27) is suppressed at least by a factor ~ 10"^"^ or 
so when compared to (ytg^h"^)^^. From figure 5 of [15] one finds (llgw/i^)_^.4 ^ 10~^^ for the linear regime in X4'^ 
preheating. Using this in (27), we see that the amplitude is hopelessly below the lower limit of observational sensitivity 
i^gwh^ ~ 10~^^ (BBO). This conclusion holds even when differences in the occupation numbers between our case and 
the Xcj)"^ model are taken into account. Assuming that both models would exhibit similar effective growth rates fj, 
but the number of oscillations would be larger by a factor AA^ in the curvaton case, the occupations numbers in the 
curvaton preheating would be enhanced by a factor exp(47r/xA7V) as compared to the Xcj)"^ case. This would enhance 
the amplitude (27) by exp(87r/iAA^) but it would still be difficult to compensate for the suppression {m„/MpY in 
(27). For example, for the case ii ~ 0.07 shown in Fig. (1) one would need A7V > 37, which is more than the number 
of oscillations taking place during the linear stage, to climb above the threshold Vlg„h? ~ 10~^*. We therefore expect 
that no directly observable gravitational waves will be produced in the linear regime of the curvaton preheating. 

However, the non-linear stages that follow the linear part of the resonance can give rise to significant production 
of gravitational waves and the dominant contribution to the total energy density of gravitational waves is typically 
generated during this stage [14, 15]. Moreover, the non-linear stages are particularly important in the curvaton model 
since the backreaction effects always become important at the final stages of preheating as discussed above. According 
to [15], the gravitational wave amplitude from these stages as observed today can be estimated as 

ng^h" ^ lQ-\RpHpf , (28) 

where i?p is the characteristic physical scale where parametric resonance is efficient and H^ is the expansion rate at 
the time when it occurs. Setting Rp ^ 1/fc, and using equations (3) and (4) we find 
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where we have defined x = m/Hp and also used the condition ga* ^ iJ* in the last step. The parameter x is 
proportional to the number of oscillation cycles and contributes a factor of few tens. The amplitude of the gravitational 
waves is therefore typically large enough to be in the detectable range, provided that the mass of the curvaton is 
not vastly smaller than the infiationary scale if*. Curvaton models with a sufficiently low infiationary scale, required 
to fulfill the frequency condition (23), might therefore produce a detectable signal of stochastic gravitational waves. 
This is an interesting outcome in light of supersymmetric models [22], for example, where the natural curvaton mass 
would be around the TeV scale. If the infiationary scale is not too much higher, these models could yield directly 
observable gravitational waves. However, in order to really confirm the qualitative outcome as well as to find out the 
exact relation between properties of the gravitational waves and the curvaton mass and infiationary scale, one would 
need to resort to a numerical study. 

DISCUSSION 

Up to now, discussions of curvaton dynamics have mostly focused on perturbative decay of the curvaton condensate. 
Here we point out that if the curvaton is about to decay directly into bosonic degrees of freedom, it is highly likely 
that there will be a non-perturbative preheating stage before the perturbative decay and reheating. The curvaton 
preheating is qualitatively similar to the conventional infiationary preheating although minor differences can arise 
due to the fact that the curvaton does not need to dominate the universe during the preheating stage. We have also 
found that the end of the curvaton preheating is always dominated by the effects of backreaction of the amplified 
light bose fields. This is different from the infiationary case, where the importance of the backreaction depends on the 
parameters. Although we have not discussed the details of backreaction, we expect to find no significant differences as 
compared to the infiationary case here either. This is because the possible differences between the curvaton preheating 
and the infiationary preheating stem from differences in the expansion of the universe which are unimportant for the 
final stages of preheating. However, to complete the analysis, one should study numerically the thermalization of 
the universe after the end of preheating. Since the curvature perturbations generated in the curvaton scenario are 
sensitive to the evolution history before the curvaton decay, these stages might potentially affect even the gaussian 
perturbations and therefore yield interesting constraints on the curvaton model itself. 

Numerics is also needed to resolve the issue of stochastic gravitational waves. Based on simple estimates we have 
postulated that models with a sufficiently low infiationary scale could produce gravitational waves with directly 
detectable frequencies and amplitudes. For example, supersymmetric models [22] with the curvaton mass around the 



TeV scale could give rise to observable gravitational waves provided that the inflationary scale is not several orders 
of magnitude higher. However, our results are still preliminary and a full numerical simulation, subject to a future 
research, is required before any firm conclusions can be drawn. 

Likewise, the possible enhancement of the non-gaussianities in curvaton preheating is best addressed by numerical 
simulations. Since curvaton preheating is very similar to the inflaton preheating, we may expect enhancements 
similar as has been found numerically by Chambers and Rajantie [13]. It therefore seems that curvaton decay by 
a parametric resonance, in addition to being a natural stage in the curvaton scenario, could also have interesting 
observational consequences. 
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